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S ' Abstract 

O ■ Due to recent important work of Zyczkowski and Sommers (J. Phys. A 36, 10115 [2003] and 

D ■ 

O , 36, 10083 [2003]), exact formulas are available (both in terms of the Hilbert-Schmidt and Bures 

metrics) for the (n 2 — l)-dimensional and ( n ( T ^~ 1 ) _ l)-dimensional volumes of the complex and 

£SJ ' real n x n density matrices. However, no comparable formulas are available for the volumes 

> 

. (and, hence, probabilities) of various separable subsets of them. We seek to clarify this situation 

O ; 

for the Hilbert-Schmidt metric for the simplest possible case of n = 4, that is, the two-qubit 

ON 
O 

O 



systems. Making use of the density matrix (p) parameterization of Bloore (J. Phys. A 9, 2059 
[1976]), we are able to reduce each of the real and complex volume problems to the calculation of 



a one-dimensional integral, the single relevant variable being a certain ratio of diagonal entries, 
v = ■ The associated integrand in each case is the product of a known (highly oscillatory 



near v = 1) jacobian and a certain unknown univariate function, which our extensive numerical 
(quasi-Monte Carlo) computations indicate is very closely proportional to an (incomplete) beta 
function B u {a, b), with a = |, 6 = \/3 in the real case, and a = ^p, b = in the complex case. 
Assuming the full applicability of these specific incomplete beta functions, we undertake separable 
volume calculations. 
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I. INTRODUCTION 

In a pair of major, skillful papers, making use of the theory of random matrices 
Sommers and Zyczkowski were able to derive explicit formulas for the volumes occupied by 
the d = (n 2 — l)-dimensional convex set of n x n (complex) density matrices (as well as the 
d = — 1 ^ TO+2 ) -dimensional convex set of real (symmetric) n x n density matrices), both in 
terms of the Hilbert-Schmidt (HS) metric [2| - inducing the flat, Euclidean geometry - 
and the Bures metric (cf. jj]). Of course, it would be of obvious considerable quantum- 
information-theoretic interest in the cases that n is a composite number, to also obtain 
HS and Bures volume formulas restricted to those states that are separable — the sum 
of product states — in terms of some factorization of n {5]. Then, by taking ratios - 
employing these Sommers-Zyczkowski results — one would obtain corresponding separability 
probabilities. (In an asymptotic regime, in which the dimension of the state space grows to 
infinity, Aubrun and Szarek recently concluded [fj] that for qubits and larger-dimensional 
particles, the proportion of the states that are separable is superexponentially small in the 
dimension of the set.) 

In particular, again for the 15-dimensional complex case, n = 4 = 2 x 2, numerical 
evidence has been adduced that the Bures volume of separable states is (quite elegantly) 
2-15(3^1) ~ 4.2136 • 10" 6 3, Table VI] and the HS volume (5V3)~ 7 w 2.73707- 10" 7 Q, eq. 
(41)]. Then, taking ratios (using the corresponding Sommers-Zyczkowski results), we have 
the derived conjectures that the Bures separability probability is 1680 ^~ 1 ) ~ 0.0733389 and 
the HS one, considerably larger, 22 ' 3 ' 7 *^ 6 13 ^ « 0.242379 {§], eq. (43), but misprinted as 5 3 
not 5 4 there]. (Szarek, Bengtsson and Zyczkowski — motivated by the numerical findings 

of aa 

- have recently formally demonstrated "that the probability to find a random 
state to be separable equals 2 times the probability to find a random boundary state to be 
separable, provided the random states are generated uniformly with respect to the Hilbert- 
Schmidt (Euclidean) distance. An analogous property holds for the set of positive-partial- 
iranspose states for an arbitrary bipartite system" [lOj (cf. [HI]). These authors also noted 
10, p. L125] that "one could try to obtain similar results for a general class of multipartite 
systems" . In this latter vein, preliminary numerical analyses of ours have given some [but 
certainly not yet conclusive] indication that for the three-qubit inseparable states, there may 
be an analogous probability ratio of 6 — rather than 2.) 
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However, the analytical derivation of (conjecturally) exact formulas for these HS and 



131 ] ) separable 



Bures (as well as other, such as the Kubo-Mori 1^] and Wigner-Yanase 
volumes has seemed quite remote — the only analytic progress to report so far being cer- 
tain exact formulas when the number of dimensions of the 15-dimensional space of 4 x 4 
density matrices has been severely curtailed (nullifying or holding constant most of the 15 



parameters) to d < 3 



HQ (cf. Q). 



Most strikingly, in this research direction, in 15, 



Fig. 11], we were able to find a highly interesting/intricate (one-dimensional) continuum 
(— oo < (3 < oo) of two-dimensional (the associated parameters being 6 1; the mean, and a^, 
the variance of the Bell-CHSH observable) HS separability probabilities, in which the golden 
ratio [17| was featured — serving to demarcate different separability regimes — among other 
items. (The associated HS volume element — 32/3(1 1 §) d 0db q da^ — is independent of b\ and 
a 2 in this three-dimensional scenario.) Further, in 14j, building upon work of Jakobczyk 



181 ]. we obtained a remarkably wide-ranging variety of exact HS separability 



and Siennicki 

(n = 4, 6) and PPT (positive partial transpose) (n = 8, 9, 10) probabilities based on two- 
dimensional sections of sets of (generalized) Bloch vectors corresponding to n x n density 
matrices. 

Nevertheless, computations for the full d = 9 and/or d = 15, n = 4 real and complex 
two-qubit scenarios are quite daunting — due to the numerous separability constraints at 
work, some being active [binding] in certain regions and in complementary regions, inactive 
[nonbinding] . "T he g eometry of the 15-dimensional set of separable states of two qubits is not 
easy to describe" [10, p. L125]. We seek to make substantial progress in these directions here, 
and, in fact, prove able to recast both these problems within one-dimensional frameworks. 



We accomplish this dimensional reduction through the use o 
of parameterization of the density matrices put forth by Bloore 



re (quite simple) form 



19 



20] some thirty years 



ago. (Of course, there are a number of other possible parametrizations 121 
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26. 2 




to estimate 



several of which we have also utilized in various studies 
volumes of separable states. Our greatest progress at this stage, in terms of increasing 
dimensionality, though, has been achieved with the Bloore parameterization — due to a 
certain computationally attractive feature of it, allowing one to decouple diagonal and non- 
diagonal parameters — which is described in sec. HT1 ) 
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A. Outline of paper 

In sec. [UJ immediately below, we describe the Bloore density matrix parameterization. 
Then, we present in sec. II III the specific one-dimensional integration formulas we have ob- 
tained for the real and complex HS separable qubit-qubit volumes using the Bloore pa- 
rameterization. The integrands in each of these cases are the product of a known jacobian 
function and a heretofore uncharacterized function. In sec. IIV| we detail the extensive numer- 
ical (quasi-Monte Carlo) procedures employed to estimate these unknown functions. Then, 
in sec. [V] we demonstrate — quite unanticipatedly — that our estimates of these functions 
over the unit interval are remarkably well-fitted (up to proportionality constants) by cer- 
tain specific incomplete beta functions. In the complex case, we can perform the indicated 
separable-volume integration exactly, but only numerically in the real case. In sec. IVI| we 
give some concluding remarks. 

In any case, it appears that further research is called for, to formally establish or appro- 
priately qualify the role of the incomplete beta function in the determination of the real and 
complex two-qubit Hilbert-Schmidt separable volumes. 



II. BLOORE PARAMETERIZATION OF DENSITY MATRICES 



The main presentation of Bloore 19( was made in terms of the 3 x 3 (n = 3) density 
matrices. It is clearly easily extendible to cases n > 3. The fundamental idea is to scale the 
off-diagonal elements (pij,i ^ j) of the density matrix in terms of the square roots of the 
diagonal entries (pa). That is, one sets (introducing the new [Bloore] variables z^), 

Pij = y/ PiiPjj Zij. (1) 

This allows the determinant of p (and analogously all its principal minors) to be expressible 
as the product (\p\ = AiA 2 ) of two factors, one (Ai = Hf =1 pu) of which is itself simply the 
product of (nonnegative) diagonal entries (pa). In the real n = 4 case under investigation 
here — we have 

A2 = (^34 — l) Z12 + 2 (Zu { Z 2A — ^23^34) + 213 (Z23 ~ Z24Z34)) Z\2 — ^23 — Z 24 ~ ^34+ (2) 
Z 14 { Z 23 ~ + Z 13 { Z 24 ~ l) + 2^23^24^34 + 2Z13Z14 (234 — ^23^24) + 1, 
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involving (only) the Zy's (i > j), where Zji = 2y [19j, eqs. (15), (17)]. Since, clearly, the 
factor A\ is positive in all nondegenerate cases (pa > 0), one can — by only analyzing A 2 

- essentially ignore the diagonal entries, and thus reduce by (n — 1) the dimensionality 
of the problem of finding nonnegativity conditions to impose on p. This is the feature we 
have sought to maximally exploit above. A fully analogous decoupling property holds in the 
complex case. 

It is, of course, necessary and sufficient for p to serve as a density matrix (that is, an 
Hermitian, nonnegative definite, trace one matrix) that all its principal minors be nonneg- 
ative [30}. The (necessary — but not sufficient) condition — quite natural in the Bloore 
parameterization — that all the principal 2x2 minors be nonnegative requires simply that 
— 1 < < 1, i ^ j. The joint conditions that all the principal minors be nonnegative are 
not as readily apparent. But for the 9-dimensional real case n = 4 — that is, Q(pij) = 

- we have been able to obtain one such set, using the Mathematica implementation of 



the cylindrical algorithm decomposition (CAD) 31[. (The set of solutions of any system of 



real algebraic equations and inequalities can be decomposed into a finite number of "cylin- 



drical" parts |32fl.) Applying it, we were able to express the conditions that an arbitrary 
9-dimensional 4x4 real density matrix p must fulfill. These took the form, 

Z12, Z13, zu G [— 1, 1], Z23 g [Z 23 , Z 23 ], z 2 4 G [Z 24 , Z 24 \, 2 34 G [Z u , Z u ], (3) 

where 



Z 23 = z 12 z 13 ± J 1 - z\ 2 J 1 - zf 3 , Z 2A = z 12 z 14 ±Jl-zf 2 Jl- zf 4 , (4) 



^34 



± Z13Z14 — Zi 2 Zi4Z 2 3 — Z 12 Z 13 Z 24: + Z 23 Z 24: ± S 



-12 



and 

S = \J—1 + ^12 + zf 3 — 1Z\ 2 Z\ 3 Z 23 + Z 23 ^J—1 + z\ 2 + zf 4 — 2Zi 2 Zi4Z 2 4 + 2^4 • (5) 

Making use of these results, we were able to confirm via exact symbolic integrations, the 
(formally demonstrated) result of Zyczkowski and Sommers [2] that the HS volume of the real 
two-qubit (n = 4) states is ~ 0.0016106. (We could also verify this through a somewhat 
[superficially, at least] different Mathematica computation, using the implicit integration 
feature first introduced in version 5.1. That is, the only integration limits employed were 
that z^ G [—1, 1], i ^ j — broader than those yielded by the CAD given by ([3]) — while the 
Boolean constraints were now imposed that the determinant of p and one [all that is needed 
to ensure nonnegativity] of its principal 3x3 minors be nonnegative.) 



A. Determinant of the Partial Transpose 



However, when we tried to combine these CAD integration limits ([3]) with the (Peres- 



Horodecki 



33, 



M 



35( n — 4) separability constraint that the determinant (A 3 = \ppr\) of 



the partial transpose of p be nonnegative 3y, Thm. 5], we exceeded the memory availabili- 
ties of our workstations. In general, the term A3 — unlike the earlier term A2 — unavoidably 
involves the diagonal entries (pa), so the dimension of the accompanying integration prob- 
lems must increase, it would seem, we initially thought — in the 9-dimensional real n = 4 
case from six to nine. 



1. Role of univariate ratio of diagonal entries 

However, we then noted that, in fact, the dimensionality of the required integrations for 
the separable volumes must only essentially be increased by one (rather than three) from 
that for the total volumes, since A3 turns out to be (aside from the necessarily nonnegative 
factor of Ax, which we can ignore) expressible solely in terms of the (six, in the real case) 
distinct Zij's and the (univariate) ratio 

v=™±. (6) 
P22P33 

(Numerical probes of ours demonstrated that v is not a local invariant of two-qubit mixed 
states, in the sense of Makhlin [371.) We, then, have 



23 I ■ 



A3 = \p PT \ = A x ( - z\^v 2 + 2z 1A (z 12 zi3 + Z24Z34) v Zt2 + su + 2z 23 (z 12 z 24 + Z13Z34) v 1/2 - z. 

(7) 

where 

S = { Z 34 ~ l) Z l 2 ~ ^ (^14^23 + ^13^24) ^34^12 — ^13 + ^14^23+ { Z 13 ~~ l) ^24 ~ ^34 — % Z X3ZX4Z 2 3Z 2 4 + 1 . 



A 3 is, thus, a quartic /biquadratic polynomial in terms of \Jv (cf. 23j, |38|). (Clearly, the 
difficulty of the two-qubit separable-volume problem under study here is strongly tied to 
the high [fourth] degree of A3 in y/v. By setting either zx4 = or z 23 = 0, the degree of A 3 



can be reduced to 2 (cf. |20|)-) In the complex case, A 3complex — which we do not explicitly 
present here — also assumes the form of a quartic polynomial in \Jv. So one must deal, in 
such a setting, with 13-dimensional integration problems rather than 15-dimensional ones. 



2. attempted seven-fold exact integration 

The problem of determining the separable volumes can, thus, be seen to hinge on (in the 
real case), a seven-fold integration involving the six (independent) z^-'s and v. However, such 
requisite integrations, allowing v to vary (or even holding v constant at various values, such 
as v = 1, thus, reducing to six-fold integrations), did not appear — rather frustratingly, we 
must admit — to be exactly/symbolically performable (using version 5.2 of Mathematica). 



III. REDUCTION TO CWE-DIMENSIONAL PROBLEMS 



Making use of the "Bloore parameterization" 19) (sec. [TT] immediately above) of density 
matrices, which allows the decoupling of diagonal entries from non-diagonal entries in certain 
relevant determinant calculations, one can show that the problem of computing the 15- 
dimensional volume (y sep / comp i ex ) of the separable two-qubit systems is reducible to a one- 
dimensional integration of the form, 

Vsep/complex 2 / J ^ aC comp i ex {]/)F ccrmp i ex {l/}(lv '. (8) 

Jo 

(We measure volume in terms of the Euclidean/Hilbert-Schmidt/Frobenius norm, and 
slightly modify our notation in 39j, to indicate that we have changed from the variable 
p used there to v = p? here. The variable v — as noted earlier ([6]) — is simply a specific 
ratio p ^ p ^ of diagonal entries pa of the corresponding 4x4 density matrices p.) 

Similarly, the 9-dimensional Hilbert-Schmidt volume of the separable real density matrices 
(those with entries restricted to the real numbers) can be expressed as 

V S ep/real = 2 Ja,C Teal {v) F real (v)dv. (9) 

Jo 

The two (highly oscilla tory near v = 1) jacobian functions {J ac rea iiy) and J ac comp i ex {y)) are 



both explicitly known 



39j, sec. III.B], that is (Fig.HJ, 



z/ 3 / 2 (12 {y{y + 2) (z/ 2 + Uv + 8) + 1) log (y/U) - 5 (5z/ + 32i/ 3 - 32u - 5)) 
Jac real {v) - 3780(z/ -l) 9 

(10) 



and (Fig. Wj 

Jar i„Jv\ = — 

3603600 (u 



1/ 



;noenn/ „ T ^ (h 1 + h 2 ), (11) 
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FIG. 1: Plot of the jacobian function Jac rea i{y), given by (fTU|) 

Jac com pi ex (v) 




FIG. 2: Plot of the jacobian function J ' ac comp i ex {v) , given by (jlip 

where 

h x = 363v 7 + 9M7u 6 + 48363// + 42875// - 42875// - 48363// - 9947// - 363; 
h2 = -140 (u 7 + 49// + Ulu 5 + 1225// 4 + 1225// 3 + 441z/ 2 + 49z/ + l) log (y/v) . 

We obtained the jacobian functions Jac rea i(u) and J (ic comp i ex {y) , given in ffTDl and ffTll . 
by transformations of, say, p 33 to the v variable (and subsequent two-fold exact integrations 
over pn and p 2 2) of the original (three-dimensional) jacobians, involving the diagonal entries, 
for the Bloore parameterizations. These original jacobians were of the form (Uf =1 pu) k with 
k = | in the real case, and k = 3, in the complex case. (Of course, by the unit trace 
condition, we must have pu = 1 — pw — P22 — P33.) 

The (only) unknowns in our two separable- volume-computation problems (([81) and (|9~1)) 
are, then, the functions F rea i(v) and i ? C omp;ei:(^)- in our preprint [39j, we reported our 
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initial numerical (quasi-Monte Carlo) procedures to estimate these two centrally important 
functions (but in terms of the variable \x = sfv). We have since continued these efforts, 
which we now detail in the following section. 



IV. ESTIMATION OF UNKNOWN UNIVARIATE FUNCTIONS 

At an advanced stage of our numerical analyses, the initial results of which had been 
reported in [39], it appeared that it might be more efficacious to employ v = p ^ p ^ as the 
principal variable rather than fi = y f^^ ■ (Thus, as previously noted, we have v = fi 2 .) 
So, our estimation (uniform sampling) procedures were originally designed in terms of 
rather than v. 

I — I 

We had, in [39], begun proceeding along two fully parallel courses, one for the 9- 
dimensional real two-qubit case and the other for the 15-dimensional complex case. We 
sought those functions freai(fi) and f complexify) , that we now see satisfy the equivalences, 

freal ( v^) = J 'real (^) ; f complex(^J~^) = ^cxmiplvxi^) : (^^) 

that would result from imposing the conditions that the expressions A\, A2 and A3 (as 
well as a principal 3x3 minor of p), along with their complex counterpart expressions, be 
simultaneously nonnegative. (The satisfaction of all these joint conditions ensures that we 
are dealing precisely with separable 4x4 density matrices.) It was evident that the relation 
/(A*) = /(~) m ust hold, so we only numerically studied the range \x G [0, 1]. Dividing this 
unit interval into 2,000 equal nonoverlapping subintervals of length each, we sought to 
estimate the /(/x)'s at the 2,001 end points of these subintervals. 

This required {p, being alternately fixed at every single one of these end points for each 
set of Zij's) numerical integrations in 6 and 12 dimensions. For this purpose, we utilized the 



Tezuka-Faure (TF) quasi-Monte Carlo procedure [40j, |4l|, we have extensively used in our 
earlier studies of separability probabilities [3, Q] (cf. [42] for an apparently more efficient 
approach to estimating the Euclidean volume of convex bodies). For each of the 2,001 
discrete, equally-spaced values of \x we employed the same set of 611,500,000 Tezuka-Faure 
six-dimensional points in the real case and, similarly, the same set of 549,500,000 twelve- 
dimensional points in the complex case. (The Tezuka-Faure points are defined over unit 
hypercubes [0, l] n , so in our computations, we transform the Bloore variables accordingly 
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and take into account the corresponding jacobians.) 



A. Close comparison with Zyczkowski-Sommers known values 

In the real case, our sample estimate of the known Hilbert-Schmidt volume of (separable 
plus nonseparable) states j^, « 0.0016106 was smaller by only a factor of 0.999996. So, 
we would expect our companion estimates of f realty), at each of the 2,001 sampled points, 
to be roughly equally precise. (Let us note that f rea i(0) = f complexity = 0- Statistical testing 
the use of confidence limits — is not appropriate in the Tezuka-Faure framework (cf. 



43J).) In the complex case, our estimate of the known 15- dimensional volume, 



851350500 



1.12925 • 10~ 6 was larger only by a factor of 1.00009. 

As instances of specific values (avoiding the necessity for 2,000 repetitions for each point), 
based on independent analyses using still larger numbers of TF-points, we obtained es- 
timates of f real (l) = F real (l) = if|||f » 114.62351, / reaZ (|) = « 74.10608, 

both based on 3,596,000,000 TF-points, and fcompiexty = F comp i ex {l) = 387.5080921 and 
fcom P iex{\) = 180.7173447, both based on 2,036,000,000 TF-points. We have the predicted 
values G complex {l) » 387.486102 and G rea i{l) « 114.6270015. (Searches using the "Plouffe's 



Inverter" website http:/ /pi.lacim.uqam.ca/eng/ did not readily yield any underlying expla- 



nation of these values or a number of transformations of them.) 

V. FITTED INCOMPLETE BETA FUNCTIONS 

Numerical computations (detailed in sec. IIVI above) - - provided us with estimates of 
Fcompiexiy) and F rea i(v) (though the sampling [quasi-Monte Carlo] procedure employed had 
been devised in terms of the variable \i = \fv and counterpart functions f C om P iex(^) and 
freai(n))- We have been able to fit these results quite well (Fig. Hj) using (concave) functions 
of the form (Fig. [3]), 



G realty) 



(^MMMM (13) 



and 



I 100000000 \ /2V6 3\ U D ^ 3 , 
G com piexty - [ jB^—^j B v ^— ^ ) . (14) 
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G(V) 




0.2 0.4 0.6 0. 



FIG. 3: The two fitted scaled incomplete beta functions G comp i ex {y) and the (lesser- valued at 
v = 1) G rea i{v) 

(Let us note that « 1.73205, ^ « 0.979796 and -| « 2.12132.) Here 5 denotes the 



(complete) beta function, and S y the incomplete beta function |44j , 



-B„(a, 6) 



w a - x (l 



w 



\b-l 



dw. 



(15) 



(The beta function itself, that is B(a,b) = BAa,b), played an important instrumental role 



in the original formulation of string theory 45|, p. 6]. For a more specific-still incomplete 



beta function role in string theory, pertaining to the symmetric group S3 and the mod- 
ular group M(2), see the review MR512916 of A. O. Barut in the MathSciNet database 
http: / / ams.rice.edu/mathscinet /| of a [somewhat obscure] paper of M. Zaganescu 



To obtain the two residual curves shown in Fig. H] — upon which we draw our central 
conclusion that F real (v) and F comp i ex (i/) are well fitted by G real {v) and G comp i ex (v), respec- 
tively — we interpolated the Tezuka-Faure points, using third-order polynomials — and 
then reparameterized the resulting curve in terms of v. 



A. separable volume and hyperarea estimations 

From the (exact) formulas of Zyczkowski and Sommers Q for the Hilbert-Schmidt vol- 
umes of the real and complex n x n density matrices to the case n = 4, we know that the 
total volume of separable and nonseparable two-qubit systems is 8513 ^ 6 0500 ~ 1.12925- 10~ 6 in 
the 15-dimensional (complex) case and ~ 0.0016106 in the 9-dimensional (real) case. 

Also, from the results of Zyczkowski and Sommers [2|, eq. (6.5)], one can readily deduce 
that the ratio of boundary (14-dimensional) hyperarea to volume of the 15-dimensional 
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FIG. 4: Our numerical (interpolated) estimates (sec. IIVP of F comp i ex (v) and F rea i(v) minus the 
values predicted by G comp i ex {u) and G rea i(v). The more strongly fluctuating curve corresponds to 
the complex case. Note the greatly reduced y-axis scale vis-a-vis that of Fig. [3J this observation 
constituting the basis for our central assertion that the i ? (i^)'s are well fitted by the £r(z/)'s. 



convex set of 4 x 4 density matrices is equal to 30a/3, and further |2j, eq. (7.9)] that the 
corresponding (lesser) ratio for the 9-dimensional convex set of real 4x4 density matrices is 
18\/3. By the subsequent results of Szarek, Bengtsson and Zyczkowski [10] — which were 
motivated by certain numerical analyses of Slater [8] — we know the analogous hyperarea- 
volume ratios for the 15- and 9-dimensional separable subsets must be simply twice as large 
(that is, 60 a/3 and 36 a/3)- 

Using the proposed incomplete beta function fits (113"]) and (jT4|) . we have attempted the 
evaluations of the two corresponding separable volumes ((jHJ) and Q), as well as separa- 
ble (lower-dimensional) hyperareas, obtaining exact results in the complex case, but only 
numerical ones for the real scenario. We succeeded in the complex case, using integration- 
by-parts, first integrating Jac comp i ex (u). (In the real case, an analogous initial integration 
of Jac rea i{y) led to a much more complicated result, now involving various hypergeometric 
functions. So, the integration by parts was stymied there.) The exact result itself in the 
complex case (for which we thank M. Trott) was very lengthy (much too so to present here), 
but we could evaluate it to any given precision. 
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1. complex case 



Using this exact formula, we were able to obtain V sep / comp i ex ~ 2.73827578 • 10 7 , and 



again applying the Zyczkowski-Sommers [2j and Szarek-Bengtsson-Zycskowski results [10], 
Psep i 'complex ~ 0.24248582 and H sep / comp i ex « 0.0000142285, all assuming the full applicabil- 
ity/validity of ffT4j) . 

We had previously hypothesized that V sep / comp i^. = (5\/3)~ 7 ~ 2.73707 • 10~ 7 [8j, eq. 
(41)] and P sep / comp iex = ^'^V 13 ^ ~ 0.242379 [sj, eq. (43), but misprinted as 5 3 not 
5 4 there]. The analysis in [8] was based on 400,000,000 quasi- Monte Carlo [Tezuka-Faure] 
points. (Those points were 15-dimensional in nature vs. the 12-dimensional ones used here.) 
Additionally, each point there was employed only once for the Peres-Horodecki separability 
test, while each point here is used in 2,000 such tests (with fx ranging over [0,1]). We had 
initially suspected that if we started checking the Peres-Horodecki criterion for successively 
larger values of fx, holding the set of Zy-'s given by a Tezuka-Faure point fixed, then if we 
reached one value of fx for which separability held, then all higher values of fx (less than 
or equal to 1) would also yield separability. But this interestingly turned out not to be 
invariably the case. So, it appeared that we needed to check the criterion 2,000 times for 
every single 6-dimensional (real) or 12-dimensional (complex) TF-point. 



2. real case 



Since, as noted, exact integration-by-parts did not seem feasible in the real case, we chose 
to expand G rea i(u) in a 75-term power series about v = 0, and performed exact integrations 
term-by-term. The overail result can be expressed as « 0.0007310253. Using 

the various results of Zyczkowski and Sommers [2J, and Szarek, Bengtsson and Zyczkowski 
[k3] detailed above, we then immediately have the estimates P sep /reai ~ 0.4538838 for the 
separability probability of the real 4x4 density matrices (markedly greater than in the 
complex case), and H sep / rea i ps 0.02279111 for the hyperarea of the bounding 8-dimensional 
hypersurface. 
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VI. CONCLUDING REMARKS 



In this study, using the Bloore parameterization of density matrices ( 19|, sec. [TTJ) , we 
have shown that incomplete beta functions (sec. E}, or clearly quite close relatives to them, 
appear to play important roles in the calculation of the Hilbert-Schmidt separable volumes 
of the 9-dimensional real and 15-dimensional complex qubit-qubit pairs. However, there 
are still apparently systematic (sine-like) — although quite small — variations (Fig. H]) of 
the estimated function from the hypothesized one G comp i ex (i>) in the complex case, so we 
suspect that we may have possibly not yet fully explained this scenario. So, to summarize, 
although we have developed here a rather compelling case for the relevance of the incomplete 
beta functions, our evidence for this is so far essentially empirical/numerical rather than 
theoretical. 

The extension to qubit-qutrit pairs (and possibly higher- dimensional composite systems, 
n > 4) of the tmwanate-function-strategy we have pursued above for the case of qubit-qubit 
pairs (n = 4), seems problematical. In the n = 4 case, the analysis is facilitated by the 
fact that it is sufficient that the determinant of the partial transpose of a density matrix be 
nonnegative for the Peres-Horodecki separability criterion to hold 36|, Thm. 5] 47]] . More 
requirements than this single one are needed in the qubit-qutrit scenario — even though 
the criterion of nonnegativity of the partial transpose is still both necessary and sufficient 
for 6x6 density matrices. (In addition to the determinant, the leading minors and/or the 
individual eigenvalues of the partial transpose of the 6x6 density matrix would need to be 
tested for nonnegativity, as well. Also the qubit-qutrit analogue of the ratio (u) of diagonal 
entries, given by (jSJ), would have to be defined, if even possible.) 

In our earlier study [20j], we had also employed the Bloore parameterization of the two- 
qubit (and qubit-qutrit) systems to study the Hilbert-Schmidt (HS) separability probabilities 
of specialized systems of less than full dimensionality. We also reported there an effort to 
determine a certain three- dimensional function (somewhat in contrast to the one-dimensional 
functions F rea i(u) and F comp i ex (u) above, but for a rather similar purpose) over the simplex 
of eigenvalues that would facilitate the calculation of the 15-dimensional volume of the 
complex two-qubit systems in terms of (monotone) metrics — such as the Bures, Kubo-Mori, 
Wigner-Yanase,. . . — other than the (non-monotone Q]) Hilbert-Schmidt one considered 
here. (The Bloore parameterization 19| did not seem immediately useful in this monotone 
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metric context, since the eigenvalues of p are not explicitly expressed (cf. 49|). Therefore, 



we had recourse in 



20(1 to the Euler-angle parameterization of density matrices of Tilma, 



Byrd and Sudarshan [23(|, in which the eigenvalues of p do, in fact, explicitly enter.) 

So, it would seem to appear, initially at least, that the particular utility of the Bloore 
parameterization in reducing the dimensionality of the problem of computing the Hilbert- 
Schmidt separable volume of qubit-qubit pairs, of which we have taken advantage in this 
study, neither extends to higher- dimensional Hilbert-Schmidt volumes (n > 4) nor to mono- 
tone metric volumes of even qubit-qubit pairs (n = 4) themselves. 
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